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Abstract 

In Combinatorial Public Projects, there is a set of projects that may be undertaken, and a set of self- 
interested players with a stake in the set of projects chosen. A public planner must choose a subset of 
these projects, subject to a resource constraint, with the goal of maximizing social welfare. Combinato- 
rial Public Projects has emerged as one of the paradigmatic problems in Algorithmic Mechanism Design, 
a field concerned with solving fundamental resource allocation problems in the presence of both selfish 
behavior and the computational constraint of polynomial-time. 

We design a polynomial-time, truthful-in-expectation, (1 — 1 /e)-approximation mechanism for wel- 
fare maximization in a fundamental variant of combinatorial public projects. Our results apply to com- 
binatorial public projects when players have valuations that are matroid rank sums (MRS), which en- 
compass most concrete examples of submodular functions studied in this context, including coverage 
functions, matroid weighted-rank functions, and convex combinations thereof. Our approximation fac- 
tor is the best possible, assuming P ^ NP. Ours is the first mechanism that achieves a constant factor 
approximation for a natural NP-hard variant of combinatorial public projects. 



'Extended abstract appears in Proceedings of the 12th ACM Conference on Electronic Commerce (EC), 2011. 
t Supported by NSF Grant CCF-0448664. 



1 Introduction 



The overarching goal of algorithmic mechanism design is to design computationally efficient algorithms 
that solve or approximate fundamental resource allocation problems in which the underlying data is a priori 
unknown to the algorithm. A problem that has received much attention in this context — albeit mostly in 
the form of negative results — is Combinatorial Public Projects (CPP). Here, there are m projects being 
considered by a public planner, n players, and a bound k <mon the number of projects that may be chosen. 
Each player i has a private valuation Vi (S) for each subset S of the projects. We consider the flexible variant 
of CPP, where a feasible solution is a set of at most k project^ The goal is to choose a feasible set of 
projects S maximizing social welfare: J2i '^i{S). The valuations are initially unknown to the public planner, 
and must be elicited from the (self-interested) players. A "mechanism" for CPP extracts this information, 
and decides on a set of projects to undertake. The mechanisms we consider can charge the players payments 
in order to incentivize truthful reporting of their valuations. Moreover, we seek mechanisms that run in 
polynomial time. 

Since CPP is highly inapproximable for general valuations — even by non-truthful algorithms — it is 
most interesting to study CPP for restricted classes of valuations. Most notable among these are submodular 
valuations, as they naturally model the pervasive notion of "diminishing marginal returns". In this paper, 
we study CPP for a fundamental and large subset of submodular valuations: Matroid Rank Sum Valuations. 
This class includes most concrete examples of submodular functions studied in this context. Most notably, 
it includes the canonical and arguably most natural example of submodularity: coverage functions. 

Combinatorial public projects and its variants are examples of welfare maximization problems. There are 
many other examples, most notable among them are combinatorial auctions, with their many variants (see 
e.g. 1 25 1). Welfare maximization problems occupy a central position in mechanism design, not only because 
of the fundamental nature of the utilitarian objective, but also due to the rich economic theory surrounding 
them. Most notably, the celebrated Vickrey-Clarke-Groves (VCG) mechanism (see e.g. [25]) is a general 
solution for all these problems, at least from an economic perspective. The VCG mechanism is truthful, in 
that it is in a player's best interest to report his true valuations regardless of the reports of the other players. 
Moreover, VCG finds the welfare maximizing solution. 

Unfortunately, however, most interesting welfare maximization problems, such as combinatorial public 
projects, are NP-hard. Therefore, implementing VCG efficiently — i.e. in polynomial time — is impossible 
unless P = NP. Moreover, as first argued in [24], most existing approximation algorithms — unlike 
exact algorithms — cannot be converted to truthful mechanisms by the imposition of a suitable payment 
scheme. This necessitates the design of carefully crafted approximation algorithms, tailored specifically 
for truthfulness. Understanding the power of these truthful approximation mechanisms is the central goal 
of algorithmic mechanism design. This research agenda was first advocated by Nisan and Ronen [23 1. 
Since then, combinatorial auctions and combinatorial public projects have emerged as the paradigmatic 
"challenge-problems" of the field, with much work in recent years establishing upper and lower-bounds on 
truthful polynomial-time mechanisms for these problems, for example: ll20l[rTl[T^[T2l[T0l l6l [l4ll27l l3ll4l 

The "holy grail" of algorithmic mechanism design is to design polynomial-time truthful approximation 
mechanisms that match the approximation guarantee of the best (non-truthful) polynomial-time approxima- 
tion algorithm. Unfortunately, several recent impossibility results have shed serious doubt on the possibility 

'This is in contrast to the exact variant, where each feasible solution consists of exactly k projects — a difference that is 
uninteresting in an approximation algorithms context, yet has major implications when incentives are in the picture. For more on 
the distinction between the two variants, we refer the reader to jTj. 
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of this goal lfT0ll27l [3ll4ir7]|. Combinatorial public projects, in particular, bore the brunt of the most bru- 
tal of these negative results ll27l 141171. Fortunately, all but one of these lower bounds apply exclusively to 
deterministic mechanisms, and none apply to randomized mechanisms for the — arguably more natural — 
flexible variant of combinatorial public projects. 

As the limitations of deterministic mechanisms became apparent, a recent research direction has focused 
on designing randomized approximation mechanisms for the fundamental problems of algorithmic mecha- 
nism design [20l IH [151 |9j |T7l|. These mechanisms are instances of the only general approach^ known for 
designing (randomized) truthful mechanisms: via maximal-in-distributional range (MIDR) algorithms [8]. 
An MIDR algorithm fixes a set of distributions over feasible solutions — the distributional range — inde- 
pendently of the valuations reported by the self-interested participants, and outputs a random sample from 
the distribution that maximizes expected (reported) welfare. The "Vickrey-Clarke-Groves (VCG)" payment 
scheme renders an MIDR algorithm truthful-in-expectation — that is, a player unaware of the coin flips of 
the mechanism maximizes his expected utility by reporting truthfully. 

Recently Dughmi, Roughgarden and Yan fTTl presented the most general framework to date for the 
design of maximal-in-distributional-range algorithms. Their approach is based on convex optimization, and 
generalizes the celebrated linear-programming based approach of Lavi and Swamy |[20l . Given a mathemat- 
ical relaxation to a welfare maximization problem, ifTTl advocates designing randomized rounding schemes 
that are convex. Given a convex rounding scheme, the problem of finding the best output of the rounding 
scheme is a convex optimization problem solvable in polynomial time, and implements an MIDR allocation 
rule. They then show how to design a convex rounding scheme for combinatorial auctions with matroid rank 
sum valuations, yielding an optimal (1 — 1/e) approximation mechanism. We elaborate on the framework 
of im in Section O 

By reducing the problem of designing a truthful mechanism to that of designing a convex rounding 
scheme, the approach of ifTTl yielded the first optimal truthful mechanism for a variant of combinatorial 
auctions with restricted valuations. It is now natural to wonder if their approach is applicable to other 
welfare maximization problems. In particular, can the convex rounding framework be used to obtain optimal 
approximation mechanisms for interesting variants of Combinatorial Public Projects? 

We answer this question in the affirmative, and elaborate on our contributions below. 

1.1 Contributions 

We design a (1 — l/e)-approximate convex rounding scheme for combinatorial public projects with matroid 
rank sum valuations. This yields a (1 — l/e)-approximate truthful-in-expectation mechanism for CPP, 
running in expected polynomial-time. This is the best approximation possible for this problem, even without 
truthfulness, unless P = NP. Therefore, ours is the first truthful mechanism for an NP-hard variant of CPP 
that matches the approximation ratio of the best non-truthful algorithm. Our results works with "black-box" 
valuations, provided that players can answer a randomized analogue of value oracles. 

To prove our results, we follow the general outline of fTTl. However, our task is more challenging: 
whereas in combinatorial auctions, randomized rounding may allocate each item independently (the ap- 
proach taken in UT]), this is not possible in CPP. We must respect the cardinality constraint of k on the set 
of chosen projects, and therefore our rounding scheme must by fiat be dependent. This presents a major 
challenge in analyzing our rounding scheme. Whereas the expected value of a submodular function on a 
product distribution (i.e. independent rounding) has been studied extensively, and is closely related to the 

^The random sampling approach used in |6|, while arguably general, does not seem applicable beyond auction settings — in 
particular, it is not applicable to combinatorial public projects. 
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now well-understood multi-linear (see e.g. ||5]|30l), analyzing the expected value of a dependent distribution 
— in particular proving it to be a concave function of underlying parameters — is a technical challenge that 
we overcome by combining techniques from combinatorics, convex analysis, and matroid theory. 

1.2 Additional Related Work 

Combinatorial Public Projects, in particular its exact variant, was first introduced by Papadimitriou, Schapira 
and Singer ll27l . They show that no deterministic truthful mechanism for exact CPP with submodular val- 
uations can guarantee better than a 0(-y/m) approximation to the optimal social welfare. The non-strategic 
version of the problem, on the other hand, is equivalent to maximizing a submodular function subject to 
a cardinality constraint, and admits a (1 — l/e)-approximation algorithm due to Nemhauser, Wolsey and 
Fisher ||2TI . and this is optimal |[28l assuming P ^ NP. 

Buchfuhrer, Schapira and Singer [4| explored approximation algorithms and truthful mechanisms for 
CPP with various classes of valuations in the submodular hierarchy. The most relevant result of [4 1 to our 
paper is a lower-bound of 0{^/m) on deterministic truthful mechanisms for the exact variant of CPP with 
coverage valuations — a class of valuations for which our randomized mechanism for flexible CPP obtains 
a (1 — 1/e) approximation. 

Most recently, Dobzinski fj] showed two lower bounds for CPP in the value oracle model: A lower 
bound of 0{^/m) on universally truthful mechanisms for flexible CPP with submodular valuations, and a 
lower bound of 0{^/m) on truthful-in-expectation mechanisms for exact CPP with submodular valuations. 
We note that the latter was the first unconditional lower bound on truthful-in-expectation mechanisms. 

2 Preliminaries 

2.1 Combinatorial Public Projects 

In Combinatorial Public Projects there is a set [m] = {!,..., m} of projects, a cardinality bound k such that 
< k < m, and a set [n] = {1, . . . , n} of players. Each player i has a valuation function Vi : 21™] — )• M+ 
that is normalized (vj(0) = 0) and monotone ivi{A) < Vi{B) whenever A C B). In this paper, we 
consider the. flexible variant of combinatorial public projects: a feasible solution is a set S* C [m] of projects 
with |5| < k. Player Ts value for outcome S is equal to Vi{S). The goal is to choose the feasible set S 
maximizing social welfare: Vi{S). 

We consider Combinatorial Public Projects where each player's valuation vi is know to lie in some set 
V of valuation functions. We abbreviate the set of instances of CPP constrained to valuations V as CPP(V). 
As first defined in [27], CPP was considered with V equal to the set of monotone submodular functions. 
In this paper, we focus on CPP with matroid-rank-sum (MRS) valuations — a large subset of monotone 
submodular functions. 

2.2 Mechanism Design Basics 

We consider direct-revelation mechanisms for combinatorial public projects. Fix m,n, and k, and let S = 
{S C [m] : 151 < A;} denote the set of feasible solutions. A mechanism comprises an allocation rule, which 
is a function from (hopefully truthfully) reported valuation functions ui, . . . , ti„ : 2["^] — )• M to a feasible 
outcome S £ S, and a payment rule, which is a function from reported valuation functions to a required 
payment from each player. We allow the allocation and payment rules to be randomized. 
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A mechanism with allocation and payment rales A and p is truthful-in-expectation if every player always 
maximizes its expected payoff by trathfully reporting its valuation function, meaning that 

nv^{A{v)) - pi{v)] > B[Vi{AU,V.^))-pM,V-i)] (1) 

for every player i, (true) valuation function Vi, (reported) valuation function v'^, and (reported) valuation 
functions v-i of the other players. The expectation in ^ is over the coin flips of the mechanism. 

The mechanisms that we design can be thought of as randomized variations on the classical VCG mecha- 
nism, as we explain next. Recall that the VCG mechanism is defined by the (generally intractable) allocation 
rule that selects the welfare-maximizing outcome with respect to the reported valuation functions, and the 
payment rule that charges each player i a bid-independent "pivot term" minus the reported welfare earned 
by other players in the selected outcome. This (deterministic) mechanism is truthful; see e.g. f2T\. 

Now let dist{S) denote the probability distributions over the feasible set S, and let V C dist{S) be a 
compact subset of them. The corresponding Maximal-In-Distributional-Range (MIDR) allocation rule is de- 
fined as follows: given reported valuation functions vi, . . . ,Vn, return an outcome that is sampled randomly 
from a distribution D* £ V that maximizes the expected welfare Es^D[J2i'^i{^)] o^^i" distributions 
D £ V. Analogous to the VCG mechanism, there is a (randomized) payment rale that can be coupled with 
this allocation rule to yield a truthful-in-expectation mechanism (see L8.I). 

2.3 Matroid Rank Sum Valuations 

We now define matroid rank sum valuations. Relevant concepts from matroid theory are reviewed in Ap- 
pendix IB. II 

Definition 2.1. A set function v : 2["^1 M is a matroid rank sum (MRS) function if there exists a family 
of matroid rank functions ui, . . . , : 2['"1 — )• M, and associated non-negative weights wi, . . . ,Wk G 
such that v{S) = ^£=i 'WeUi{S) for all S C [m]. 

We do not assume any particular representation of MRS functions, and require only oracle access to their 
(expected) values on certain distributions (see Section l2!4l ). MRS valuations include most concrete examples 
of monotone submodular functions that appear in the literature — this includes coverage function^ matroid 
weighted-rank function^ and all convex combinations thereof. Moreover, as shown in [28], 1 — 1/e is the 
best approximation possible for CPP with coverage valuations — and hence also for MRS valuations — 
in polynomial time, even ignoring strategic considerations. That being said, we note that some interesting 
submodular functions — such as some budget additive function^ — are not in the matroid rank sum family. 

2.4 Lotteries and Oracles 

A value oracle for a valuation v : 2["^] — > R takes as input a set S" C [m], and returns v{S). We define an 
analogous oracle that takes in a description of a simple lottery over sets S C [m], and outputs the expectation 
of V over this lottery. 

''a coverage function / on ground set [m] designates some set 3^, and m subsets y4i, Am, C 3^, such that /( 5) — \UiizsAi\. 
We note that 3^ may be an infinite, yet measurable, space. Coverage functions are arguably the canonical example of a submodular 
function. 

This is a generalization of matroid rank functions, where weights are placed on elements of the matroid. It is true, though not 
immediately obvious, that a matroid weighted-rank function can be expressed as a weighted combination of matroid (unweighted) 
rank functions — see e.g. 11161 . 

^ A set function / on ground set [m] is budgeted additive if there exists a constant B > Q (the budget) such that f{S) = 
min(B,E,es/({i}))- 
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Let k G [m], let R C [m], and let x G [0, l]'" be a vector such that J2j ^ 1- We interpret x as 
a probability distribution over [m] U {*}, where * represents not choosing a project. Specifically, project 
j G [m] is chosen with probability xj, and * is chosen with probability 1 — Xj. We define a distribution 
D^{x) over 2'™!, and call this distribution the k-bounded lottery with marginals x and promise R. We 
sample S ~ Dj}{x) as follows: Let ji, . . . , jf^ be independent draws from x, and let S = RU {ji , ■ ■ ■ ,jk}\ 
{*}. Essentially, this lottery commits to choosing projects R, and adds an additional k projects chosen 
randomly with replacement from distribution x. When = 0, as will be the case through most of this 
paper, we omit mention of the promised set. We can now define a randomized analogue of a value oracle 
that returns the expected value of a bounded-lottery. 

Definition 2.2. A bounded-lottery-value oracle /or set function v : 2[™] — )■ M takes as input a vector 
X £ [0, 1]*" with ^ ^> <^ bound k G [m], and a set R C [m], and outputs ^q^j~,r/^\[v{S)]. 

In our model for CPP, we assume that a player with valuation function Vi can answer bounded-lottery- 
value oracle queries for Vi. A bounded-lottery-value oracle is a generalization of value oracles. Nevertheless, 
it is the case that a bounded-lottery-value oracle can be implemented using a value oracle for some succinctly 
represented examples of MRS valuations, such as explicit coverage functions (In similar fashion to ifTT] 
Appendix A]). 

More generally we note that bounded-lottery-value oracles can be approximated arbitrarily well, with 
high probability, using value oracles; this is done by random sampling, and we omit the technical details. 
Unfortunately, we are not able to reconcile the incurred sampling errors — small as they may be — with 
the requirement that our mechanism be exactly truthful. We suspect that relaxing our solution concept to 
approximate truthfulness - also known as e-truthfulness - would remove this difficulty, and allow us to relax 
our oracle model to the more traditional value oracles. 

2.5 Convex Rounding 

In this section, we review convex rounding, a framework for the design of truthful mechanisms introduced 
by Dughmi, Roughgarden and Yan [.17 J . We present the main definitions and lemmas as they pertain to 
combinatorial public projects. For a more thorough and general treatment of convex rounding, we refer the 
reader to [ 17 Section 3]. 

We consider the standard integer programming formulation of CPR There is a variable Xj G {0, 1} 
for each project j G [m], and the goal is to set at most k of the variables to 1 so that the welfare v{x) = 
Vi{{j : Xj = 1}) is maximized. We relax this integer program in the obvious way to the polytope V = 

|x G M'" : Y.-Xj <k,xt o}. We postulate a rounding scheme r that maps points of V to the feasible 
solutions S = {S CI [m] : |5| < k} of CPP. We allow r to be randomized, so that r{x) is a distribution over 
S for each x £ V. 

Traditionally, approximation algorithms optimize an objective v{x) — often a simple extension of v to 
V — over the set V of fractional solutions, and then round the optimal fractional point x* to a solution r{x*) 
in the original feasible set S. Many of the best approximation algorithms for various problems are based 
on this relax-solve-round framework. Unfortunately, however, this approach is almost always incompatible 
with the design of truthful mechanisms, due to the fact that the rounding step is often unpredictable. Truthful 
mechanism design, on the other hand, is intimately tied to exact optimization, as evidenced by the fact that 
the vast majority truthful mechanisms for multi-parameter problems are based on the VCG paradigm (see 
Section |22]l. 
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In an effort to reconcile the techniques of approximation algorithms and truthful mechanism design, 
Dughmi, Roughgarden and Yan proposed optimizing directly on the output of the rounding scheme, rather 
than on its input. This defines an optimization problem induced by relaxation V and rounding scheme r. 
Stated for CPP with the relaxation as described above, the problem is as follows. 



They consider a simple allocation rule, which we state for CPP in Algorithm [U that solves dlj) optimally. 
They observe that this allocation rule is maximal-in-distributional-range. 

Algorithm 1 MIDR Allocation Rule for CPP 
Parameter: n,m,k 

Parameter: (Randomized) rounding scheme r 
Input: Valuation functions {vi}^^^ 
Output: A set C [m] with \S\ < k 
1: Let X* be an optimal solution to Q 

2: Let S ~ r{x*) 



Lemma 2.3 ( ifTTl '). Algorithm\l\is an MIDR allocation rule. 

For a < 1, we say that the rounding scheme r for CPP(V) is a-approximate if, whenever x is an integer 
point of V corresponding to a set 5" G S, and G V for each i, we have that £2-^^(3,) • fj(T)] > 
a Vi{S). In other words, rounding does not degrade the quality of an integer solution by more than a. 
Given the definition of Algorithm [B it is easy to conclude the following lemma. 

Lemma 2.4 ([17]). If r is an a-approximate rounding scheme for CPP(V), then Algorithm \J\ is an a- 
approximation algorithm for CPP(V). 

For reasons outlined in lITTl . implementing Algorithm [T] efficiently is impossible for most rounding 
schemes r in the literature. To get around this difficulty, they advocate designing rounding schemes that 
render dD a convex optimization problem. 

Definition 2.5. Consider a randomized rounding scheme r : "P — )• dist{S). We say r is a convex rounding 
scheme /or CPP(V) if whenever vi G V for all i, the objective E5^f,(j.)[^,- Vi{S)] is a concave function of 

X. 

Lemma 2.6. When r is a convex rounding scheme for CPP(V) , ^ is a convex optimization problem for 
each instance of CPP(V). 

Under additional technical conditions, discussed in the context of combinatorial public projects in Ap- 
pendix El convex program (|2]i can be solved efficiently (e.g., using the ellipsoid method). This reduces the 
design of a polynomial-time a-approximate MIDR algorithm to designing a polynomial-time a-approximate 
convex rounding scheme. 

Summarizing, Lemmas 12.31 12.41 and 12.61 give the following informal theorem. 

Theorem 2.7 (Informal). If there exists an a-approximate convex rounding scheme for CPP(V), then there 
exists a truthful-in-expectation, polynomial-time, a-approximate mechanism for CPP(V). 




for j = 1, . . . 



m. 



(2) 
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3 The Mechanism 



In this section, we prove the main result. 

Theorem 3.1. Thereisa {1 — 1/ e)-approximate, truthful-in-expectation mechanism for combinatorial public 
projects with matroid rank sum valuations in the bounded-lottery-value oracle model, running in expected 
poly(n, m) time. 

We structure the proof of Theorem 13. li as follows. We define the k-bounded-lottery rounding scheme, 
which we denote by r^, in Section |3?T] We prove that is (1 — 1/e) -approximate (Lemma 1331 ). and convex 
(Lemma 13.21 ). Lemmas 12.31 12.41 and 13.31 taken together, imply that Algorithm [H when instantiated with 
r = r/c, is a (1 — l/e)-approximate MIDR allocation rule. Lemma [l!2] reduces implementing this allocation 
rule to solving a convex program. 

In Appendix |A1 we handle the technical and numerical issues related to solving convex programs. First, 
we prove that our instantiation of Algorithm [T]can be implemented in expected polynomial-time using the 
ellipsoid method under a simplifying assumption on the numerical conditioning of our convex program 
(Lemma IA.2I ). Then we show in Section IA.3I that the previous assumption can be removed by slightly 
modifying our algorithm. 

Finally, we prove that truth-telling VCG payments can be computed efficiently in Lemma IB. 41 Taken 
together, these lemmas complete the proof of Theorem 13.11 

3.1 The /c-Bounded- Lottery Rounding Scheme 

We devise a rounding scheme r^. that we term the k-bounded-lottery rounding scheme. Given a feasible 
solution X to linear program (|2), we let distribution rk{x) be the /c-bounded-lottery with marginals x/k (and 
promise 0), as defined in Section 1241 We make this more explicit in Algorithm |2] 

Algorithm 2 The A;-Bounded-Lottery Rounding Scheme 

Input: Fractional solution x G with Xj < k, and < < 1 for all j. 

Output: S C [m] with |5| < A: 

1: For each j £ [m] designate the interval Ij = J2j'<j ^j'l i ^j'<j ^i'] of length ^ 

2: Draw pi, . . . ,pk independently and uniformly from [Ci, 1] 

3: Let S = {j e [m] : {pi, ...,pk}nlj^ 0} 



The A;-bounded-lottery rounding scheme is (1 — 1/e) approximate and convex. We prove the approxima- 
tion lemma below. As for convexity, we present a simplified proof for the special case of coverage valuations 
in Section |3^ and present the proof for MRS valuations in Section [331 

Lemma 3.2. The k-bounded-lottery rounding scheme is convex for CPP with MRS valuations. 

Lemma 3.3. The k-bounded-lottery rounding scheme is (1 — l/e)-approximate when valuations are sub- 
modular 

Proof. Fix n,m,k and {vi}^^^. Let S C [m] be a feasible solution to CPP — i.e. \S\ < k. Let be 
the vector with 1 in indices corresponding to S, and otherwise. Let T ~ ^^(15). We will first show that 
each element of j G S" is included in T with probability at least 1 — 1/e. Observe that T is the union of 
k independent draws from a distribution on [m] U {*}, where each time the probability of j € S is 1/k. 
Therefore, the probability that j is included in T is 1 — (1 — 1/k)'' > 1 — 1/e. 
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Submodularity now implies that E[fj(r)] > (1 — 1/e) • Vi{S) for each player i — this was proved in 
many contexts: see for example |[T9l Lemma 2.2], and the earher related result in lITSl Proposition 2.3]. This 
completes the proof. □ 



3.2 Warmup: Convexity for Coverage Valuations 

In this section, we prove a special case of Lemma lB^ for coverage valuations. Recall that a coverage function 
/ on ground set [m] designates some set 3^, and m subsets Ai, . . . , A„i ^ 3^, such that f{S) = \ Uj^s 

Fix n,m,k and {viY^^^. Assume that, for each player i, the valuation function Vi : 2^"^^ — )• M is a 
coverage function. We let v{S) = J2i ^i('S') be the welfare of a solution S to CPP. It is an easy observation 
that the sum of coverage functions is also a coverage function. Therefore v{S) is a coverage function. We 
let 3^ be a set, and Ai, . . . , Am ^ y, such that v{S) = \ Uj^s Aj\. While our proof extends easily to the 
case where y is an arbitrary measure space, we assume in this section that 3^ is a finite set for simplicity. 

Let V denote the polytope of fractional solutions to CPP as given in We now show that Eg^j-^ix) [^(5*)] 
is a concave function of x for a; G V, completing the proof of Lemma [3^ for the special case of coverage 
valuations. Take an arbitrary x £ V, and let S ~ ^^(x) be a random variable. Using linearity of expecta- 
tions, we can rewrite the expected welfare as follows. 

E[viS)] = E[| U,^sA,\] = Y,Pr[£ G U.esAj] 

Since the sum of concave functions is concave, showing that Pr G Uj^sAj] is concave in x for each ^ G 3^ 
suffices to complete the proof. For ^ G 3^, let Ti = {j G [m] : i £ Aj} he the set of projects that "cover" 
£. Let pi, . . . ,pk and /i, . . . , be as in Algorithm|2] Note that are disjoint sub-intervals of [0, 1], 

and \Ij\ = We can rewrite the probability of covering i as follows. 

Pr[i G UjesA,] = Fr[S n / 0] 

= Pr[{pi, ... ,pk} nu jerjj ^ 9] 
= 1 - Pr[{pi, ...,pk}n Uj^Tjj = 0] 

k 

t=i 
k 

= i_jj(i_meT,/,i) 

t=i 

The final form is simply the composition of the concave function g{y) = 1 — (1 — y/k)^ with the affine 
function y — )• X^^gj-^ xj. It is well known that composing a concave function with an affine function yields 
another concave function (see e.g. IH). This completes the proof. 

3.3 Convexity for Matroid Rank Sum Valuations 

In this section, we will prove Lemma [3l2] in its full generality. First, we recall the discrete hessian matrix, 
as defined in Iil7il . 
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Definition 3.4 ( II17I ). Let v : 2^"^! ^be a set function. For S C [m], we define the discrete Hessian 
matrix V.s£W^^"^ofvatS as follows: 

nUi, j) = vis U {i, j}) - v{S U {i}) - vis U {j}) + viS) (3) 

forij G [m]. 

It was shown in ifTTl that the discrete hessian matrices are negative semi-definite for matroid rank sum 
functions. 

Claim 3.5 (Inl). Ifv : 2H ^ m+ 

is a matroid rank sum function, then Tig is negative semi-definite for 

each S C [m]. 

We now return to Lemma |3]2] Fix n and m. For each cardinality bound k £ [m], let Vj. denote the 
polytope of fractional solutions to CPP as given in For a set of MRS valuations vi, . . . , we observe 
that the social welfare viS) = Y17=i ^iiS) is — by the (obvious) fact that the sum of MRS functions is an 
MRS function — also an MRS function. Therefore, we will prove Lemma 13.21 by showing that, for each 
k G [m] and MRS function v : 2[™1 — )■ M, the following function of x G Vk is concave in x. 

Glix) = E HS)] 

= viS)Pv[rkix) = S] ^"^^ 

SC[m] 

We use techniques from combinatorics to write Pr[rkix) = 5] in a fonii that will be easier to work 
with. For T C [m], we use xt as short-hand for XljeT ^i' ^^'^ ^ short-hand for [m] \ T. 

Claim 3.6. For each k G [m], x € Vt, cind S C [m] 

Pr[r,(x) = 5] = -ll^l (5) 

Proof. It is easy to see that Pr[rjt(x) = S] is equal to: 



k 

RCS 



Pr[nix) CS]- Pr[\/ Tkix) C S \ {j}] (6) 

Using the inclusion-exclusion principle, we can rewrite (|6]l as follows: 

Pr[rfc(x)C5]- -ll^l"'Pr[rfc(x) C5\r] (7) 

Letting = 5 \ T in ©, we get 

Pr[nix) C 5] - ^ Pr[rkix) C R] (8) 

RCS 

We can easily simplify ^ to conclude that 

Pr[rkix) = S]=Y -l'^'"'^' Pr[rfc(x) C R] (9) 

RCS 
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Next, we observe that the expression Pr[rfc(x) C R] can be expressed as a simple closed form in x. 
Let pi, . . . ,pk and /i, ...,/„ be as in Algorithm |2] The event rfc(x) C R occurs exactly when none of 
pi, . . . ,Pk land in the intervals corresponding to projects R. Recalling that the interval Ij of project j has 
length Xj/k, we get that the probability of any particular pt falling in ^j^-jilj is exactly x^/k. Therefore, 
by the independence of the variables pi, . . . we get that 



Pr[rfc(x) CR]=i^l-^j (10) 

Combining ^ and ( fTOb completes the proof. □ 

Building on Claim [3^ we now express the Hessian matrix of as a non-negative weighted sum of 

. We note that when x G P^, it is eaf 

therefore ([TTI i is well-defined. 

Claim 3.7. For each k G [m], x G Vk, and v : 2['"1 — > M, we have 



discrete Hessian matrices of v. We note that when x G Vk, it is easy to verify that • x G Vk-2, and 



5C[m' 



ni (11) 



Proof. Fix j G [m], possibly with i = j. We work with as defined in Equation and plug in 
expression (|5]l. 

G^.(-)= E -(5)--iI^IE-1"''(i-t)' 

scH Rcs 

Differentiating with respect to Xj and Xj gives: 



5C[m] RQS\{i,j} 



We group the terms by projecting 5 onto [m] \{i,j}, and then we simplify the resulting expression. 



d^Glix) ^ k-1 Y_^\s\ Y_i\R\ (i 



- U {i}) - ^(5 U {j}) + U {i,i})) 



■SC[m]\{»j} «C5 

=^ E (12) 

SC[m] ACS 

The second equality follows from the fact that v{S) — v{S U {i}) — v{S U {j}) + v{S U {i, j}) = when 
S includes either of i and j. The last equality follows by definition of Tig. 

k 

([T2I ). we conclude that 

52G^(x) A; - 1 



Invoking Claim 13.61 with k' = k — 2 and x' = ■ x, and plugging the resulting expression into into 



dxidxj k 



Epi 



5C[m] 



rk^2 — ; x ] = S 



□ 
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Claims (13.51) and (13.71 ) establish that, when v is MRS and k G [m], \/'^G'^{x) is a non-negative weighted 
sum of negative semi-definite matrices for each x £ Vk- ^ non-negative weighted sum of negative semi- 
definite matrices is negative semi-definite. Therefore, the Hessian matrix of is negative semi-definite at 
each X G Vk, and we conclude that is a concave function on Vk- This completes the proof of Lemma 
EH 
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A Solving The Convex Program 



In this section, we overcome some technical difficulties related to the solvability of convex programs. We 
follow the general outline of iflTl Appendix B], modifying the proofs throughout in order to handle the 
additional technical difficulties specific to CPP. We show in Section lA.ll that, in the bounded-lottery-value 
oracle model, the four conditions for "solvability" of convex programs, as stated in Fact IB.3I are easily 
satisfied for convex program (|2]i when r = r^. However, an additional challenge remains: "solving" a 
convex program — as in Definition IB .21 — returns an approximately optimal solution. Indeed the optimal 
solution of a convex program may be irrational in general, so this is unavoidable. 

We show how to overcome this difficulty if we settle for polynomial runtime in expectation. While 
the optimal solution x* of ^ cannot be computed explicitly, the random variable (x* ) can be sampled 
in expected polynomial-time. The key idea is the following: sampling the random variable rj^{x*) rarely 
requires precise knowledge of x*. Depending on the coin flips of r^, we decide how accurately we need 
to solve convex program (|2]i in order compute rk{x*). Roughly speaking, we show that the probability of 
requiring a (1 — e)-approximation falls exponentially in -. As a result, we can sample rk{x*) in expected 
polynomial-time. We implement this plan in Section IA.2I under the simplifying assumption that convex 
program is well-conditioned — i.e. is "sufficiently concave" everywhere. In Section IA31 we show how 
to remove that assumption by slightly modifying our algorithm. 

A.l Approximating the Convex Program 

Claim A.l. There is an algorithm for Combinatorial Public Projects with MRS valuations in the bounded- 
lottery-value oracle model that takes as input an instance of the problem and an approximation parameter 
e > 0, runs in poly(n, m, log(l/e)) time, and returns a (1 — e)-approximate solution to convex program ^ 
when r = r^. 

It suffices to show that the four conditions of Fact IB. 31 are satisfied in our setting. The first three are 
immediate from elementary combinatorial optimization (see for example ||29l ). It remains to show that the 
first-order oracle, as defined in Fact IB. 31 can be implemented in polynomial-time in the bounded-lottery- 
value oracle model. We let f{x) denote the objective function of convex program (|2]) when r = r^. This 
objective can, by definition, be written as follows. 



where Vi is the valuation function of player i and G^' is as defined in (01). By definition, G'"^{x) is the 
outcome of querying the bounded-lottery- value oracle of Vi with bound k and marginals x/k. Therefore, we 
can evaluate f{x) using n bounded-lottery-value queries, one for each player. It remains to show that we 
can also evaluate the (multi-vaiiate) derivative v/(x) of /(x). Using definition (01) and Claim [321 we take 
the partial derivative of G^' with respect to Xj and simplify the resulting expression. 
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x] =S 



SC\m\ 



rk-1 



k-1 



S 



(13) 



The second equality follows by grouping the temis of the summation by the projection of S onto [m] \{j}. 
The third equality follows from the observation that v{S U {j}) — v{S) = when S includes j. The fourth 
equality follows by a simple re-arrangement and application of Claim [331 

Inspect the final form ([T3] ) in light of the definition of bounded-lottery- value oracles (Definition 12.21 ) and 
the definition of (Section [3TT] ). Notice that the first term is the expected value of Vi over the {k — 1)- 
bounded-lottery with marginals and promise {j}. The second term is the expected value of Vi over the 

same lottery without the promise. Therefore, we can evaluate (x) using two queries to the bounded- 
lottery-value oracle of player i. This completes the proof of Claim IaIT] 



A.2 The Well-Conditioned Case 

In this section, we make the following simplifying assumption: The objective function f{x) of convex 
program Q with r = r^, when restricted to any line in the feasible set V, has a second derivative of 

v-(\m\) 

magnitude at least A = ^p^M^.m) everywhere, where the polynomial in the denominator may be arbitrary. 
This is equivalent to requiring that every eigenvalue of the Hessian matrix of f{x) has magnitude at least A 
when evaluated at any point in V. Under this assumption, we prove Lemma ED 

Lemma A.2. Assume the magnitude of the second derivative of f{x) is at least A = ^piji"y^l^l^^"* everywhere. 
AlgorithmUl instantiated with r = r^, can be simulated in time polynomial in n and m in expectation. 

Let X* be the optimal solution to convex program with r = r^. Algorithm [T]outputs a set of projects 
distributed as rfc(a;*). The A;-bounded-lottery rounding scheme, as described in Algorithm |2j requires mak- 
ing k independent decisions: for £ S {1, ... , k}, we draw uniformly from [0, 1] and decide which interval 
Ij, if any, pi falls into. In other words, we find the minimum index (if any) such that J2j<j^ x*/k > p£. 
Fix i. For most realizations of pi, we can calculate j£ using only coarse estimates xj to x*. Assume we have 
an estimation oracle for x* that, on input 6, returns a 6-estimate x of x*: Specifically, xj — x* < 6 for each 
j G [m]. If pi falls outside the "uncertainty zones" of x, such as when \p£ — Ylj'<j xj' /k\ > 6m/k for each 
j G [m], it is easy to see that we can correctly determine ji by using x in lieu of x. The total measure of 
the uncertainty zones of x is at most 2m? 6, therefore p£ lands outside the uncertainty zones with probability 
at least 1 — 2ir?5. The following claim shows that if the estimation oracle for x* can be implemented in 
time polynomial in log(l/(5), then we can simulate the /c -bounded-lottery rounding procedure in expected 
polynomial-time. 
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Claim A. 3. Let x* be the optimal solution of convex program Q with r = r^. Assume access to a subroutine 
B{5) that returns a 5-estimate of x* in poly(n, m, log(l/(5)) time. AlgorithmU} instantiated with r = r^, 
can be simulated in expected poly(n, m) time. 

Proof. Fix i G {1, . . . , fc}. Draw pi G [0, 1] uniformly at random as in the A;-bounded-lottery rounding 
scheme in Algorithm |2] We will show how to find, in expected poly(n, m) time, the minimum index ji (if 
any) such that Y.j<j, ^j/^ > Pi- 

The algorithm proceeds as follows: Start with S = 6o = Let x = B{6). While \pi—Ylj'<j ^j' /^l ^ 
6m /k for some j G [m] (i.e. p£ may fall inside an "uncertainty zone") do the following: let 5 = 6/2, 
x = B{6) and repeat. After the loop terminates, we have a sufficiently accurate estimate of x* to calculate 
je- 

It is easy to see that the above procedure is a faithful simulation of Algorithm Q on x*. It remains 
to bound its expected running time. Let 6t = 2^+1^^ denote the value of 6 at iteration t. By our initial 
assumption, iteration t takes poly(n, m, log{l/ 6t)) = poly(n, m, log(2*+^m^)) = poly(n, m, t) time. The 
probability this procedure does not terminate after t iterations is at most 2m? 6t = 1/2*. Taken together, 
these two facts and a simple geometric summation imply that the expected runtime is polynomial in n and 
m. □ 

It remains to show that the estimation oracle B{6) can be implemented in poly(n, m, log(l/5)) time. 
At first blush, one may expect that the ellipsoid method can be used in the usual manner here. However, 
there is one complication: we require an estimate x that is close to x* in solution space rather than in terms 
of objective value. Using our assumption on the curvature of /(x), we will reduce finding a 5-estimate of 
X* to finding an 1 — e{6) approximate solution to convex program dU with r = r^. The dependence of e on 
6 will be such that e > poly(<5)/2P°iy("'™), thereby we can invoke Claim IaTT] to deduce that B{6) can be 
implemented in poly(n, m, log(l/5)) time. 

Let e = e{6) = 2J2^li[m]) - Plugging in the definition of A, we deduce that e > (^2y2P°iy(">'^), which is 
the desired dependence. It remains to show that if x is (1 — e)-approximate solution to Q, then x is also a 
(5-estimate of x*. 

Using the fact that f{x) is concave, and moreover its second derivative has magnitude at least A, it a 
simple exercise to bound distance of any point x from the optimal point x* in terms of its sub-optimality 
f{x*) - f{x), as follows: 

f{x*)- f{x)>^\\x-x*\\\ (14) 
Assume 5; is a (1 — e)-approximate solution to ^ with r = r^. Equation (fT4l ) implies that 

where the last inequality follows from the fact that ?;j([m]) is an upper-bound on the optimal value /(x*). 
Therefore, ||x — x*|| < 5, as needed. This completes the proof of Lemma IA.2I 



A.3 Guaranteeing Good Conditioning 

In this section, we propose a modification of the /e-bounded-lottery rounding scheme r^. We will argue 
that satisfies all the properties of established so far, with one exception: the approximation guarantee 
of Lemma 13. 3 1 is reduced to 1 — 1/e — 2^^™". Then we will show that satisfies the curvature assumption 
of Lemma IA.2I demonstrating that said assumption may be removed. Therefore Algorithm [T] instantiated 
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with r = for combinatorial public projects with MRS valuations in the bounded-lottery-value oracle 
model, is (1 — 1/e — 2~^™") approximate and can be implemented in expected poly(n, m) time. Finally, 
we show in Remark rA.4l how to recover the 2^^™" term to get a clean 1 — 1/e approximation ratio, as 
claimed in Theorem 13. II 

Let ^ = 2~2"™. We define in Algorithm |3] Intuitively, first chooses a tentative set S C [m\ of 
projects using r^. Then it cancels its choice with small probability ^. Finally, with probability /3 it chooses 
a random project j* G [m] and lets S = {j*}. /3 is defined as the fraction of projects included in the original 
tentative choice of S. The motivation behind this seemingly bizarre definition of is purely technical: as 
we will see, it can be thought of as adding "concave noise" to r^. 

Algorithm 3 Modified A;-bounded-lottery Rounding Scheme 

Input: Fractional solution x G M™ with Y^- Xj < k, and < < 1 for all j. 

Output: Feasible solution S C [m] with \S\ < k 

1: Lets' = rfc(x) 

2: Let /? = ^ 

3: Draw qi G [0, 1] uniformly 
4: if qi G [0, /x] then 
5: Let S = 

6: Draw 92 £ [0, 1] uniformly 
7: if q2 G [0, /3] then 

8: Choose project j* G [m] uniformly at random. 

9: Let S = {j*} 
10: end if 
11: end if 



We can write the expected welfare E^^^+^^^ • Vi{S)] as follows. 



E 

5'-rfc(z) 



Using linearity of expectations and the fact that (3 is independent of the choice of j* to simplify the expres- 
sion, we get that E^_^+(^) • Vi{S)] is equal to 



(1 - ^) E 



+ /iE[ 



EiLiEr=i^i({j}) 



m 



Observe that chooses a project j with probability 1 — (1 — Xj/k^. Therefore, the expectation of /3 is 



E, i-(i-x,7fc)* 



. This gives: 



E 

S~r+(x) 



(1 - //) E 



(15) 



It is clear that the expected welfare when using r = is within 1 — ^ = 1 — 2 of the expected 
welfare when using r = in the instantiation of Algorithm [T] Using Lemma [331 we conclude that is 
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a (1 — 1/e — 2~^"™)-approximate rounding scheme. Moreover, using Lemma |3^ as well as the fact that 
1 — (1 — Xj/k)^ is a concave function, we conclude that is a convex rounding scheme. Therefore, this 
establishes the analogues of Lemmas 13.31 and \i?2\ for r^. It is elementary to verify that our proof of Lemma 
IA.2l can be adapted to as well. 

It remains to show that is "sufficiently concave". This would establish that the conditioning assump- 
tion of Section IAT2I is unnecessary for r^. We will show that expression ([T5] ) is a concave function with 

V v-{\ra]) 

curvature of magnitude at least A = ^^i^znm everywhere. Since the curvature of concave functions is 
always non-positive, and moreover the curvature of the sum of two functions is the sum of their curvatures, 
it suffices to show that the second term of the sum ([15] ) has curvature of magnitude at least A. We note that 
the curvature of (l — (1 — xj/k)^) is at least over x G [0, 1]™. Therefore, the curvature of the 
second term of (031 ) is at least 

as needed. 

Remark A.4. In this section, we sacrificed 2"^"™ in the approximation ratio in order to guarantee expected 
polynomial runtime of our algorithm even when convex program ^ is not well-conditioned. This loss can be 
recovered to get a clean 1 — 1/e approximation as follows. Given our (1 — 1/e — 2~^^"^) -approximate MIDR 
algorithm A, construct the following algorithm A': Given an instance of combinatorial public projects, 
A' runs A on the instance with probability 1 — 62^^"*", and with the remaining probability solves the 
instance optimally in exponential time 0(2^"*"). It was shown in [15] that a random composition of MIDR 
mechanisms is MIDR, therefore A' is MIDR. The expected runtime of A' is bounded by the expected runtime 
of A plus 62"^""* • 0(2^"'") = 0(1). Finally, the expected approximation of A' is the weighted average of 
the approximation ratio of A and the optimal approximation ratio 1, and is at least (1 — e2^^"™)(l — 1/e — 

2-2nm) _^ g2-2"™ > 1 - 1/e. 

B Additional Preliminaries 
B.l Matroid Theory 

In this section, we review some basics of matroid theory. For a more comprehensive reference, we refer the 
reader to f26l. 

A matroid M is a pair (Af , X), where X is a. finite ground set, and X is a non-empty family of subsets of 
X satisfying the following two properties. (1) Downward closure: If S belongs to X, then so do all subsets 
of S. (2) The Exchange Property: Whenever T, S £ I with |r| < IS"], there is some x G 5 \ T such that 
T U {x} € X. Elements of X are often referred to as the independent sets of the matroid. Subsets of X that 
are not in X are often called dependent. 

We associate with matroid M a set function rankM : 2'^ — )• N, known as the rank function ofM, defined 
as follows: rankM{A) = m.axs£i\Sr\A\. Equivalently, the rank of set ^4 in matroid M is the maximum size 
of an independent set contained in A. A set function / on a ground set is a matroid rank function if there 
exists a matroid M on the same ground set such that / = rankM- Matroid rank functions are monotone 
{f{S) < f{T) when 5 C T), normalized (/(0) = 0), and submodular (/(5) + /(r) > f{SnT) + f{SuT) 
for all S and T). 
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B.2 Convex Optimization 

In this section, we distill some basics of convex optimization. For more details, see lUl. 

Definition B.l. A maximization problem is given by a set li of instances {V, c), where V is a subset of 
some euclidean space, c : "P — )• M, and the goal is to maximize c{x) over x & V. We say H is a convex 
maximization problem if for every ("P, c) € 11, "P is a compact convex set, and c : V ^ M is concave. If 
c : P — )• for every instance of H, we say U is non-negative. 

Definition B.2. We say a non-negative maximization problem 11 is i?-solvable in polynomial time if there 
is an algorithm that takes as input the representation of an instance I = (P,c) G 11 — where we use 
\I\ to denote the number of bits in the representation — and an approximation parameter e, and in time 

]io\y{\I\,\og{l/e)) outputs X such that c{x) > (1 — e) max^g-p c(y). 

Fact B.3. Consider a non-negative convex maximization problem 11. If the following are satisfied, then 11 is 
R-solvable in polynomial time using the ellipsoid method. We let X = {V,c) denote an instance o/Il, and 
let m denote the dimension of the ambient euclidean space. 

1. Polynomial Dimension: m is polynomial in 

2. Starting ellipsoid: There is an algorithm that computes, in time poly(|X|), a point c £ W^, a matrix 
A G and a number V G M such that the following hold. We use E{c, A) to denote the ellipsoid 
given by center c and linear transformation A. 

(a) E{c,A) D P 

(b) V < volumeiV) 

volume{E{c,A)) ^ 2Poly(|2^l) 

3. Separation oracle for V: There is an algorithm that takes takes input I and x £ W^, and in time 
poly(|X|, |x|) where \x\ denotes the size of the representation of x, outputs "yes" ifx G V, otherwise 
outputs h G M™ such that x < h^yfor every y £ V. 

4. First order oracle for c: There is an algorithm that takes input I and x G W\ and in time poly \x\) 
outputs c{x) G M and \jc{x) G M™". 

B.3 Computing Payments 

Lemma B.4. Let A be an MIDR allocation rule for combinatorial public projects, and let vi, . . . ,Vn be 
input valuations. Assume black-box access to A, and value oracle access to {vi}^^^. We can compute, with 
poly(n) over-head in runtime, payments pi, . . . ,pn such that 'E[pi] equals the VCG payment of player ifor 
MIDR allocation rule A on input vi, . . . ,Vn- 

We note that an essentially identical lemma was proved in lITTl . Nevertheless, we include a proof for 
completeness. 

Proof. Without loss of generality, it suffices to show how to compute pi. Let : 2'™! — M be the valuation 
evaluating to at each bundle. Recall (see e.g. 1221 ') that the VCG payment of player 1 is equal to 

E 

T~A(0,V2,...,v„) 



.1=2 



E 



.1=2 



(16) 
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Let S be a sample from A{vi, . . . , Vn), and let T be a sample from .4(0, V2, ■ ■ ■ , Vn)- Let pi = 
Yll=2 ^i(^) ~ Z]r=2 ^j('S')- Using linearity of expectations, it is easy to see that the expectation of pi is 
equal to the expression in (fT6l ). This completes the proof. □ 

We note that the mechanism resulting from Lemma lR!4l is individually rational in expectation, and each 
payment is non-negative in expectation. We leave open the question of whether it is possible to enforce 
individual rationality and non-negative payments for oiu" mechanism ex-post. 
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